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Abstract 

The quantum statistical parton distributions approach proposed more 
than one decade ago is revisited by considering a larger set of recent and 
accurate Deep Inelastic Scattering experimental results. It enables us to 
improve the description of the data by means of a new determination of the 
parton distributions. This global next-to-leading order QCD analysis leads 
to a good description of several structure functions, involving unpolarized 
parton distributions and helicity distributions, in a broad range of x and 
and in terms of a rather small number of free parameters. There are several 
challenging issues and in particular the confirmation of a large positive gluon 
helicity distribution. The predictions of this theoretical approach will be 
tested for single-jet production and charge asymmetry in production 
in pp and pp collisions up to LHC energies, using recent data and also for 
forthcoming experimental results. 
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1 Introduction 


Deep Inelastic Scattering (DIS) of leptons and nucleons is indeed our main 
source of information to study the internal nucleon structure in terms of 
parton distributions. Several years ago a new set of parton distribution func¬ 
tions (PDF) was constructed in the framework of a statistical approach of 
the nucleon [T]. For quarks (antiquarks), the building blocks are the helicity 
dependent distributions g^(a;) (g^(a;)). This allows to describe simultane¬ 
ously the unpolarized distributions q{x) = q~^{x) + q~{x) and the helicity 
distributions Aq{x) = q'^{x) — q~{x) (similarly for antiquarks). At the initial 
energy scale Q^, these distributions are given by the sum of two terms, a 
quasi Fermi-Dirac function and a helicity independent diffractive contribu¬ 
tion. The flavor asymmetry for the light sea, i.e. d{x) > u{x), observed in 
the data is built in. This is simply understood in terms of the Pauli exclu¬ 
sion principle, based on the fact that the proton contains two up-quarks and 
only one down-quark. We predict that d{x)/u{x) must remain above one for 
all X values and this is a real challenge for our approach, in particular in 
the large x region which is under experimental investigation at the moment. 
The flattening out of the ratio d[x)/u{x) in the high x region, predicted 
by the statistical approach, is another interesting challenge worth mention¬ 
ing. The chiral properties of QCD lead to strong relations between q{x) and 
q{x). For example, it is found that the well established result Au{x) > 0 
implies Au{x) > 0 and similarly Ad{x) < 0 leads to Ad{x) < 0. This ear¬ 
lier prediction was conhrmed by recent polarized DIS data and it was also 
demonstrated that the magnitude predicted by the statistical approach is 
compatible with recent BNL-RHIC data on production (see Section 5). 
In addition we found the approximate equality of the flavor asymmetries, 
namely d{x) — u{x) ~ Au{x) — Ad{x). Concerning the gluon, the unpolar¬ 
ized distribution G{x, Ql) is given in terms of a quasi Bose-Einstein function, 
with only one free parameter, and for simplicity, we were assuming zero gluon 
polarization, i.e. AG{x,Qq) = 0, at the initial energy scale. As we will see 
below, the new analysis of a larger set of recent accurate DIS data, has 
forced us to give up this assumption. It leads to the confirmation of a large 
positive gluon helicity distribution, giving a significant contribution to the 
proton spin, a major point which was emphasized in a recent letter [2]. In 
our previous analysis all unpolarized and helicity light quark distributions 
were depending upon eight free parameters, which were determined in 2002 
(see Ref. m), from a next-to-leading (NLO) £t of a small set of accurate DIS 
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data. Concerning the strange quarks and antiquarks distributions, the sta¬ 
tistical approach was applied using slightly different expressions, with four 
additional parameters [3]. Since the first determination of the free parame¬ 
ters, new tests against experimental (unpolarized and polarized) data turned 
out to be very satisfactory, in particular in hadronic reactions, as reported 
in Refs. 

It is crucial to note that the quantum-statistical approach differs from the 
usual global parton fitting methodology for the following reasons: 

i) It incorporates physical principles to reduce the number of free parameters 
which have a physical interpretation 

ii) It has very specific predictions, so far confirmed by the data 

iii) It is an attempt to reach a more physical picture on our knowledge of 
the nucleon structure, the ultimate goal would be to solve the problem of 
confinement 

iv) Treating simultaneously unpolarized distributions and helicity distribu¬ 
tions, a unique siuation in the literature, has the advantage to give access to 
a vast set of experimental data, in particular up to LHC energies 

The paper is organized as follows. In Section 2, we review the main points of 
our approach and we describe our method to determine the free parameters 
of the PDF with the set of experimental data we have used. In Section 3, 
we exhibit all the unpolarized and helicity distributions we have obtained. 
In Section 4, we show the results obtained for the unpolarized DIS structure 
functions Ff’“*( 0 ;, in a wide kinematic range, compared with the world 

data. We also consider e^p neutral and charged current reactions, and v{v)p 
charged current reactions. This will be completed by our analysis of polar¬ 
ized DIS experiments, like double helicity asymmetries on a proton and on 
a neutron target. In Section 5, we present predictions for cross sections and 
helicity asymmetries in hadronic collisions, in particular inclusive single-jet 
production and W production in pp and pp collisions, up to LHC energies. 
We give our final remarks and conclusions in the last section. 
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2 Basic review on the statistical parton dis¬ 
tributions 


Let us now recall the main features of the statistical approach for building 
up the PDF, as opposed to the standard polynomial type parameterizations 
of the PDF, based on Regge theory at low x and on counting rules at large 
X. The fermion distributions are given by the sum of two terms, a quasi 
Fermi-Dirac function and a helicity independent diffractive contribution: 


xq’^{x,Ql) = 


A 


+ 


AgX^’> 


exp[(x — Xq)/x] + 1 exp(a;/a;) + 1 ’ 


( 1 ) 


xq'^{x,Ql) = 




Oq 


X 


+ 


AaX^ 


exp[(a: + XQq)/x\ + 1 exp(a;/a;) + 1 ’ 


( 2 ) 


at the input energy scale Qq = IGeV^. We note that the diffractive term is 
absent in the quark helicity distribution Ag and in the quark valence contri¬ 
bution q — q. 

In Eqs. (1T|^ the multiplicative factors Xq^ and in the numerators 

of the non-diffractive parts of the g’s and g’s distributions, imply a modih- 
cation of the quantum statistical form, we were led to propose in order to 
agree with experimental data. The presence of these multiplicative factors 
was justihed in our earlier attempt to generate the transverse momentum de¬ 
pendence (TMD) [7], which was revisited recently |S] . The parameter x plays 
the role of a universal temperature and are the two thermodynamical po¬ 
tentials of the quark g, with helicity h = ±. They represent the fundamental 
characteristics of the model. Notice the change of sign of the potentials and 
helicity for the antiquarks 0. 

For a given flavor g the corresponding quark and antiquark distributions in¬ 
volve the free parameters, X^, Aq, Aq, Aq, bq, bq and bq, whose number is 
reduced to seven by the valence sum rule, f (q(x) — q(x))dx = Nq, where 
Nq = 2,1,0 for u, d, s, respectively. 

For the light quarks g = u,d, the total number of free parameters is 
reduced to eight by taking, as in Ref. d], A^ = Ad, ^ = Ad, A^ = Ad, bu = 
bd, bu = bd and b^ = bd- For the strange quark and antiquark distributions. 


^ At variance with statistical mechanics where the distributions are expressed in terms 
of the energy, here one uses x which is clearly the natural variable entering in all the sum 
rules of the parton model. 
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the simple choice made in Ref. [T] was improved in Ref. [3], but here they 
are expressed in terms of seven free parameters. 

For the gluons we consider the black-body inspired expression 


xG{x,Ql) 


Agx’’^ 

exp{x/x) — 1 ’ 


(3) 


a quasi Bose-Einstein function, with bo being the only free parameter, since 
Ag is determined by the momentum sum rule. In our earlier works IB El , we 
were assuming that, at the input energy scale, the polarized gluon, distribu¬ 
tion vanishes, so 

xAG{x,Ql) = 0. (4) 

However as a result of the present analysis of a much larger set of very ac¬ 
curate unpolarized and polarized DIS data, we must give up this simplifying 
assumption. We are now taking 


xAG{x,Ql) 


Agx^(^ 1 

(1 -h cgx'^g'^ exp{x/x — 1) 


(5) 


It is clear that we don’t have a serious justification of the functional form of 
AG{x, Ql). However the above expression shows that it is strongly related to 
G{x, Ql) and therefore constructed by means of a Bose-Einstein distribution 
with zero potential. Actually since AG{x, Ql) = P{x)G{x, Ql) a simpler 
expression would be P{x) = Ax^, but the additional term x'^^ in the denom¬ 
inator is needed in order to get a reasonable fit of the polarized DIS data. 
To insure that positivity is satisfied we must have |T’(a;)| < 1 (see section 
3). However for quarks and antiquarks positivity is automatically fullfied by 
construction. 

To summarize the new determination of all PDF’s involves a total of twenty 
one free parameters: in addition to the temperature x and the exponent bG 
of the gluon distribution, we have eight free parameters for the light quarks 
{u,d), seven free parameters for the strange quarks and four free parame¬ 
ters for the gluon helicity distribution. These parameters will be determined 
from a next-to-leading order (NLO) QCD £t of a large set of accurate DIS 
data, unpolarized and polarized structure functions, as we will discuss in the 
following section. 
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3 Unpolarized and polarized parton distribu¬ 
tions 

In order to determine these parameters we have performed a global NLO 
QCD htting procedure using only DIS data, because it is well known that 
the consideration of semi-inclusive DIS data involves uncertainties related to 
fragmentation functions. For unpolarized DIS we have considered Q"^) 

from NMC, E665, HI, ZEUS, neutral and charged current cross sections 
from HERA and charged current neutrino and anti-neutrino cross sections 
from CCER, NuTeV and CHORUS, which allow to extract xF^^{x, Q'^). We 
present in Table 1 the details of the number of points and corresponding 
for each experiment, with a total of 1773 data points for a total of 2288. 
For polarized DIS we have considered gl’^’"'{x, Q'^) from HERMES, E155, 
SMC, EMC, E143, E154, JLab and COMPASS. We present in Table 2 the 
details of the number of points and corresponding for each experiment, 
with a total of 269 data points for a total x^ of 319. 

The PDF QCD evolution was done in the MS scheme using the HOPPET 
program [37], the minimization of the x^ was performed with the CERN 
MINUIT program |38|. For unpolarized and polarized data we work in the 
General Mass Variable Flavour Number Scheme (GM-VFNS) [iQl |4T] and 
for the heavy quark we have taken mc=1.275 GeV. For the strong running 
coupling as{Q‘^) we took as{Ql) = 0.32 and we hnd as{Ml) = 0.119 ±0.001. 
In our calculations as{M^) is not a free parameter, it comes from the the 
evolution equations with an initial value ^^(Qo) ~ 0.32. For polarized data 
we have not introduced additional constraints coming from the hyperon decay 
constants, but we predict at = 2GeV^, F ± D = 1.23 ± 0.03 and 3F - 
D = 0.57 ± 0.02, to be compared respectively with the experimental values 
1.269 ± 0.003 and 0.586 ± 0.031 [I2]|l3]. Goncerning nuclear corrections 
for data involving nuclear targets, in the case of polarized data for gf, we 
have taken into account the correction due to D-wave. For unpolarized data, 
there is an effect for F 2 only in the high x region [39] that we have not 
considered. However for neutrino DIS data on iron target, the ratio of proton 
over neutron is included. For the kinematic cuts we have used the values 
given by the experiments and we have restricted the data to > IGeV^ 
and X > 10“^. We are aware that target mass corrections can be included 
to improve the data description [39], but in this simple approach they were 
not considered. Also we have taken only leading-twist effects. The error 
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bands were calculated using the standard Hessian matrix method, following 
the prescription described in Ref. |10] and we have used the standard choice 
of tolerance =1. 
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process 


^data 

xVd-o.f. 

X^ 2002 

da (z/ p) CCFR [9] 

271 

172 

1.57 

*410 

da {u p) NuTeV [TO] 

206 

177 

1.16 

*390 

da {y p) CHORUS [11] 

78 

64 

1.22 

*176 

da {p p) CCFR [9] 

191 

163 

1.17 

*318 

da {p p) NuTeV [TO] 

153 

125 

1.22 

*266 

E665 [12] 

24 

11 

2.18 

17 

ZEUS [13] 

26 

17 

1.53 

26 

Ff HI [Hug 

105 

70 

1.5 

*292 

NMC [IE] 

12 

14 

0.85 

13 

Fi NMC [IE] 

230 

155 

1.48 

274 

F^/F^ NMC [H] 

259 

205 

1.48 

198 

F^ - F^ NMC [H ng 

17 

9 

1.88 

49 

xF^^ CHORUS [H] 

65 

47 

1.38 

*89 

xF^^ NuTeV [IE] 

68 

49 

1.38 

100 

Charged current e^p HERA [20] 

33 

32 

1.03 

*75 

Charged current e~p HERA [20] 

18 

31 

0.58 

*24 

Neutral current e'^p HERA [20] 

343 

285 

1.20 

*520 

Neutral current e~p HERA [20] 

185 

138 

1.34 

*383 

Fl hi [21] 

5 

9 

0.56 

*37 

Total 

2288 

1773 




Table 1: Detailed for the cross sections and the unpolarized structure 
functions. In the column 2002 the values marked with an asterisk were 
not htted in 2002. 






process 


data 

xVd.o.f. 

2002 

gl HERMES [22] 

34 

34 

1 

*36 

gi E155 [23| 

6 

8 

0.75 

8 

gl SMC [211 

25 

12 

2.08 

35 

gl EMC [261 [2l| 

9 

10 

0.9 

8 

gl E143 [28] 

33 

28 

1.18 

34 

gl COMPASS [23 

16 

9 

1.77 

*23 

gl SMC [211 

4 

7 

0.57 

8 

gl E155 [23] 

10 

11 

0.91 

12 

gl E154 [30] 

5 

11 

0.45 

6 

gl E143 [2H] 

41 

27 

1.52 

43 

gl Jlab [31] 

3 

3 

1 

*1 

gf HERMES [221 |3l] 

43 

36 

1.19 

*70 

gf COMPASS [35] 

12 

10 

1.2 

*30 

gf E155 [36] 

21 

23 

0.91 

*43 

gf E143 [2E] 

34 

28 

1.21 

43 

gf SMC [21] 

23 

12 

1.92 

*37 

Total 

319 

269 




Table 2: Detailed for the polarized structure functions In 

the column 2002 the values marked with an asterisk were not htted in 
2002. 


The new determination of the pdf| leads, for the light quarks (g = u,d), 
to the following parameters: 

Ag = 1.943 ± 0.005, bq = 0.471 ± 0.001 Ag = 8.915 ± 0.050, 

bq = 1.301 ± 0.004, Ag = 0.147 ± 0.003, bq = 0.0431 ± 0.003 (6) 

^ To compute the unpolarized distributions and the helicity distributions, a Fortran 
program is available upon request. 
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and four potentials 


X+^ = 0.475 ± 0.001, Xq" = 0.307 ± 0.001, 

X+ = 0.245 ± 0.001, X- = 0.309 ± 0.001. (7) 

Concerning the strange quarks we have the following parameters: 

A, = 28.508 ± 0.005, 6, = 0.370 ± 0.002, A, = 0.0026 ± 0.0002, 

b, = 0.201 ± 0.003, i, = 13.689 ± 0.050, 6, = 9.065 ± 0.020, (8) 

and two potentials 

Xo+ = 0.011 ± 0.001, Xg- = 0.015 ± 0.001. (9) 

Finally in the gluon sector, we obtain the following parameters: 

Ag = 36.778 ± 0.085, be = 1.020 ± 0.0014, = 26.887 ± 0.050, 

be = 0.163 ± 0.005, cg = 0.006 ± 0.0005, dc = -6.072 ± 0.350. (10) 

In addition the new universal temperature is a; = 0.090 ± 0.002. 

By comparing with the results of 2002 |1], we have observed a remarkable 
stability of some important parameters, the light quarks potentials X^ and 
X^, whose numerical values are almost unchanged B The new temperature 
is slightly lower. As a result the main features of the new light quark and 
antiquark distributions are only scarcely modihed. However it is instructive 
to note that in Tables 1 and 2, one can judge the improvement obtained in 
this new version compared to the old 2002 version. The last column gives 
the obtained with the old parameters and the new data set. 

We display in FigJT]the different unpolarized parton distributions xf{x,Q^) 
(/ = u, d, s, c, u, d, s and G) versus x, after NLO QCD evolution at = 
lOGeV^ , with the corresponding error bands. Similarly the different quark 
and antiquark helicity distributions xAf{x, Q^) (/ = u, d, s, u, d and s) ver¬ 
sus X, after NLO QCD evolution at = lOGeV^ , with the corresponding 
error bands are shown in Fig. [2J 

Our determination of the gluon helicity distribution deserves a special dis¬ 
cussion. We display in Fig. [3]Top the gluon helicity distribution versus x at 
the initial scale Qq = IGeV^ and = lOGeV^. At the initial scale it is 

^Note the interesting relation Xq^ ~ already found in Ref. [T]. 
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sharply peaked around x = 0.4, but this feature lessens after some QCD 
evolution. We note that P{x) introduced above, has the following expres¬ 
sion, P{x) = -|- 0.006), which is such that 0 < P{x) < 1 

for 0 < a; < 1, so positivity is satished and the gluon helicity distribution 
remains positive. As already mentioned the term x’^^ plays an important 
role. It has a strong effect on the quality of the fit of Q^), since the 

increases substantially when do decreases. Its value also affects the shape 
of the gluon helicity distribution, which becomes larger towards the smaller 
X-values, for smaller do- We display AG{x, Q‘^)/G{x, in Fig. ^^ottom for 
two values and some data points [HI 05], which suggest that the gluon 
helicity distribution is positive indeed. According to the constraints of the 
counting rules this ratio should go to 1 when x = 1, but we observe that this 
is not the case here, since for example at the initial scale P{x = 1) = 0.726. 
In some other parameterizations in the current literature, this ratio goes to 
zero, since the large x behavior of xAG(x) is (1 — x)^, with (3 » 3 [HI 07] . 
Clearly one needs a better knowledge of AG(x, Q^)/G{x, Q"^) for x > 0.2. 
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Figure 1: The different unpolarized parton distributions xf{x, (/ = 
M, d, s, c, u, d, s and G) versus x, after NLO QCD evolution at = lOGeV^, 
with the corresponding error bands. The charm distributions xc{x, = 
xc{x, Q‘^) are generated by QCD evolution. 
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10 ® 10 ® 10 ^ 1 


X 

Figure 2: The different quark and antiquark helicity distributions xAf{x, Q"^) 
(/ = u, d, s, u, d and s) versus x, after NLO QCD evolution at = lOGeV^, 
with the corresponding error bands. The charm helicity distributions gener¬ 
ated by QCD evolution are essentially zero 
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0.15 


= IGeV" 
= lOGeV^ 



- = 2 GeV^ 

= 10 GeV^ 

- 0.1 L-^^—' . i _ - ^^^ — . 

10 10 ^ 1 
X 

Figure 3: Top : The gluon helicity distribution xAG{x, Q^) versus x, for 
= 1 GeV^ (dashed curve) and after NLO QCD evolution for = 10 GeV^ 
(solid curve), with the corresponding error band . 

Bottom : AG(x, (5^)/G(x, Q^) versus x, for = 2 GeV^ (solid curve) and 
= 10 GeV^ (dashed curve). The data are from HERMES [H] and COM¬ 
PASS @5] 
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Deep inelastic scattering 

4.1 Unpolarized DIS experiments 

First we present some selected experimental tests for the unpolarized PDF 
by considering /iiV and eN DIS, for which several experiments have yielded a 
large number of data points on the structure functions Q^), N stands 

for either a proton or a deuterium target. We have used hxed target mea¬ 
surements which probe a rather limited kinematic region in and x and 
also HERA data which cover a very large range and probe the very low 
X region, dominated by a fast rising behavior, consistent with our diffractive 
term (See Eq. ([T])). 

For illustration of the quality of our £t and, as an example, we show in Fig. 0] 
and Fig. [5l our results for on different hxed proton targets, to¬ 

gether with HI and ZEUS data . We note that the analysis of the scaling 
violations leads to a gluon distribution xG{x,Q^), in fairly good agreement 
with our simple parametrization (See Eq. (]3|)). 

Another rather interesting physical quantity is the neutron structure 
function and in particular the ratio F^ / F 2 {x, Q‘^) which provides strong con- 
traints on the PDF of the nucleon. For example the behavior of this ratio at 
large x is directly related to the ratio of the dto u quarks in the limit x —)■ 1, 
a long-standing problem for the proton structure. We show the results of two 
experiments, NMC in Fig. |6l which is very accurate and covers a reasonable 
range up to x = 0.7 and CL AS in Fig. [TJ which covers a smaller 
range up to larger x values, both are fairly well described by the statistical 
approach. Several comments are in order. In the small x region this ratio, for 
both cases, tends to 1 because the structure functions are dominated by sea 
quarks driven by our universal diffractive term. In the high x region dom¬ 
inated by valence quarks, the NMC data suggest that this ratio goes to a 
value of the order of 0.4 for x near 1, which corresponds to the value 0.16 for 
d{x)/u{x) when x —)■ 1, as found in the statistical approach [S]. The CLAS 
data at large x cover the resonance region of the cross section and an impor¬ 
tant question is whether Bloom-Gilman duality holds as well for the neutron 
as it does for the proton. We notice that the predictions of the statistical 
approach suggest an approximate validity of this duality, except for some 
low values. A better precision and the extension of this experiment with 
the 12GeV Jefferson Lab will certainly provide even stronger constraints on 
PDFs up to X ~ 0.8. 
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Figure 4: Ff {x, Q^) as a function of x for fixed {Q'^) and data from HERMES 
gS], E665 gg, NMC gT], EMC gS], HI glgg, ZEUS g3], BCDMS gU]. 
The function c{i) = 0.6(16 —i), i = 1 corresponds to (Q^) = 1.25GeV^ . The 
curves are the results of the statistical approach. 
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Figure 5: -F|’(a;, as a function of x for fixed high (Q^) and data from HI 
[m US], ZEUS [IS]. The function c{i) = 0.6(14 — i), i = 1 corresponds to 
(Q^) = 400GeV^ . The curves are the results of the statistical approach. 
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Figure 6: Q"^) as a function of x for fixed (Q^) and data from NMC 

[51]. The function c{i) = 0.6(20 —i), i = 1 corresponds to (Q^) = 1.125GeV^. 
The curves are the results of the statistical approach. 
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Figure 7: Q^) as a function of x* ~ a; (for the exact definition of x* 

see Ref. [52]) for fixed (Q^) and data from CLAS BoNus [52]. The function 
c{i) = 0.7(10 — i), i = 1 corresponds to (Q^) = 1.05GeV^. The curves are 
the results of the statistical approach 
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We now turn to the inclusive neutral and charged current cross sec¬ 
tions which, in addition to Ff, give access to other structure functions. 

The neutral current DIS processes have been measured at HERA in a kine¬ 
matic region where both the 7 and the Z exchanges must be considered. The 
cross sections for neutral current can be written, at lowest order, as 


dxdQ'^ 


27ra‘^ 


Y+F2{x,Q^)tY.xF:,{x,Q^) 


y"FL{x,Q^) , 


( 11 ) 


where 


F2{x, Q^)=F^ ( x , Q^)-VeXz{Q^)F^^{x, Q^)^al+vl)xliQ^)F^{x, Q^) , ( 12 ) 


xF^ix, Q^) = -aeXz{Q‘^)xF^^{x, Q^) + 2aeVexliQ^)xFi{x, Q‘^) . (13) 

The structure function Fl{x, Q^) is sizeable only at high y and we will come 
back to it later. The other structure functions introduced above, have the 
following expressions in terms of the parton distributions 0 


p 7 pjZ pZ 

-^2 5^2 , ^2 

xF^^,xF^ 


{x,Q^) 

(x,Q 2 ) 


Y1 K’ «/ + ^/] Q^) + Q^)), 

/ 

[2e/a/, 2afVf] {xqf{x, Q^) - xqf{x, Q^)) . (14) 

/ 


Here the kinematic variables are y = Q^/xs, Y± = 1±(1—p)^, y/s = 2^Ef,Ep, 
Ef> and Ep are the electron (positron) and proton beam energies respectively. 
Morever, Vi and Oj are the vector and axial-vector weak coupling constants 
for the lepton e and the quark /, respectively, and ej is the charge. The 
function Xz(Q^) is given by 


Xz(Q^) 


1 

sin^ 26*w 


(15) 


where div is the weak mixing angle and Mz is the Z-boson mass. The reduced 
cross sections are defined as 




Q^x d'^a^fz 
Ki_27rQ!^ dxdQ"^ 


(16) 


Our predictions are compared with HERA data in Fig. [HI as a function of x, 
in a broad range of values and the agreement is good. 

For low the contribution of the longitudinal structure function Fl(x, Q^). 
to the cross section at HERA is only sizeable at x smaller than approximately 

For simplicity we write them at LO, but the calculations were done by including the 
NLO corrections [55] . 
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Figure 8: Top : Comparison of the data on the reduced neutral cross section 
Q^), in e~p collisions as a function of x and for different values, 
with the results of the statistical approach. Data are from HERA |20j . 
Bottom : Same for e'^p collisions. 
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10 ^ and in this domain the glnon density dominates over the sea qnark den¬ 
sity. More precisely, it was shown that |5l] 


Fl{x,Q^) 


as{Q^) 

71 


+2Sj 


[i 

[3 A y y 
Jx y y y 


(17) 


Before HERA was shnt down, a dedicated run period, with reduced proton 
beam energy, was approved, allowing HI to collect new results on Fl. We 
show on Fig. |9]the expectations of the statistical approach compared to the 
new data, whose precision is reasonable. The trend and the magnitude of 
the prediction are in fair agreement with the data, so this is another test 
of the good predictive power of our theoretical framework. One can also 



Figure 9: The longitudinal proton structure function Fl{x,Q^) averaged in 
X at given values. Data are from ZEUS [SS] and HI [SS] and the curve is 
the result of the statistical approach. 

test the behavior of the interference term between the photon and the Z 
exchanges, which can be isolated in neutral current e^p collisions at high Q^. 
We have to a good approximation, if sea quarks are ignored, xF^^{x, Q'^) = 
+ dv){x, Q^) and comparison between data and prediction is shown in 

Fig. [TOl 
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Figure 10: The interference term extracted in e^p collisions at HERA. 
Data are from ZEUS |57| and HI [58], the curve is the prediction of the 
statistical approach. 


The charged current DIS processes have been also measured accurately 
at HERA in an extented kinematic region. It has a serious impact on the 
determination of the unpolarized parton distributions by allowing a flavor 
separation because they involve only the exchange. The cross sections 
are expressed, at lowest order, in terms of three structure functions as follows 


^ '^Born 

dxdQ"^ 


Gj. M‘y 
4ir (Q^ + Ml)^ 

^Y_xF^(x,Q")\ 


[YyF^{x,Q-^) 


5 


l^FFi 


X, 


(18) 


and the reduced cross sections are dehned as 


a 


(x,Q2) = 


Gl 




[dvr (g2 + M^)2 


-1 


dxdQ'^ 


(19) 


At leading order for e p —?• with a longitudinally polarized beam 


F^ix^Q"^) = x[m(x, g^) + c(x,g^) + d(x, g^) + s(x, g^)] 
xF^'^{x,Q‘^) = x[m(x, g^) + c(x, g^) - d(x, g^) - s(x, g^)], ( 20 ) 
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and for e'^p —)■ v^X 

F 2 ''{x,Q‘^) = x[d{x,Q‘^) + s{x,Q‘^) + u{x,Q‘^) + c{x,Q‘^)] 
xF^'^{x,Q‘^) = x[d{x,Q‘^) + s{x,Q'^) - u{x,Q'^) - c{x,Q‘^)] . (21) 


At NLO in QCD is non zero, but it gives negligible contribution, except 
at y values close to 1. Our predictions for at NLO are compared 

with HI and ZEUS data in Fig. [TTl as a function of x in a broad range of 
values. 

The differential inclusive neutrino and antineutrino cross sections have the 
following standard expressions 


dxdy 


GlMpE, g2) 




2E„ 


y 


My - y 




x,Q^ 


( 22 ) 


y is the fraction of total leptonic energy transfered to the hadronic system 
and Ey is the incident neutrino energy. F 2 and F 3 are given by Eq. fl20l) for 
up and Eq. for up, and Fi is related to F 2 by 


2xFi = 


1 + Ax^M^IQ^ 

- 1 , T-. - M 


(23) 


where R = cii/cTT, the ratio of the longitudinal to transverse cross sections 
of the W-boson production. Our results at NLO compared with the CCER 
and NuTeV data are shown in Fig. [T2l As expected, for fixed x, the y depen¬ 
dence is rather flat for neutrino and has the characteristic (1 — y)"^ behavior 
for antineutrino. 
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Figure 11: Left : Comparison of the data on the reduced charged cross section 
{x,Q‘^), in e~p collisions as a function of x and for different values, 
with the results of the statistical approach. Data are from HERA |2Uj . 

Right : Same for collisions. 
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Figure 12: Comparison of the data on the differential cross sections v{v)N 
for Ey = 150 GeV, as a function of y and for different x values, with the 
results of the statistical approach. Data are from CCFR [5] and NuTeV [TU] . 
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4.2 Polarized DIS experiments 

The spin-dependent structure function gi{x,Q^) has the well-known NLO 
QCD expression [5^ 


9i{x,Q^) 




as{Q^) 

271 


6C,) 


as{Q^) 

271 


AG 


^1 


(24) 


Aq{x, Q^), Aq{x, Q^) and AG{x, Q^) are quark, antiquark and gluon helicity 
distributions in the nucleon. 6C{x)q^G are the NLO spin-dependent Wil¬ 
son coefficient functions and the symbol ® denotes the usual convolution in 
Bjorken x space, n/ is the number of active flavors for light quarks. 

We recall that according to the results shown in Section 3, we have obtained 
a good flavor separation of these helicity distributions: for all x and val¬ 
ues, Au > 0 is the largest one. Ad < 0 is smaller in magnitude, Au > 0 and 
Ad < 0 are approximately opposite and As < 0, As < 0 are much smaller. 


We now turn to the important issue concerning the asymmetries Q^), 

measured in polarized DIS. We recall the definition of the asymmetry Ai{x, Q"^), 
namely 


Ai{x,Q‘^) 


[9i{.x, - 7^(a;, Q‘^)g2{,x, Q^)] 2a:[l + R{x, Q^)] 

F2(a;,Q2) [l + j^{x,Q^)] 


(25) 


where 9 i, 2 {x, Q"^) are the polarized structure functions, 7 ^( 0 ;, Q^) = 4:X^Mp/Q‘^ 
and R{x, Q"^) is the ratio between the longitudinal and transverse photoab¬ 
sorption cross sections. When x ^ 1 for = 4 GeV^, R is the order of 
0.30 or less and 7 ^(x, Q^) is close to 1, so if the u quark dominates, we have 
Ai ~ 0.6Au{x)/u{x), so it is unlikely to find Ai —)■ 1 , as required by the 
counting rules prescription, which we don’t impose. We display in Fig. [13] 
the world data on A{{x, Q^) {Top) and A^{x, Q^) {Bottom), with the results 
of the statistical approach at = 4GeV^, up to x = 1. Indeed we hnd that 

Ar < 1. 

Finally one important outcome of this new analysis of the polarized DIS 
data in the framework of the statistical approach, is the confirmation of a 
large positive gluon helicity distribution, which gives a significant contribu¬ 
tion to the proton spin [ 2 ]. 
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Figure 13: Top : Comparison of the world data on with the result 

of the statistical approach at = 4 GeV^, including the corresponding error 
band. 

Bottom : Comparison of the world data on Q‘^) with the result of the 

statistical approach at = 4 GeV^, including the corresponding error band. 
Data are from [22] - [33] . 






















































5 Hadronic collisions 


A precise determination of parton distributions allows us to use them as 
input information to predict strong interaction processes, for additional tests 
of pertubative QCD and also for the search of new physics. Here we shall 
test our statistical parton distributions for the description of two inclusive 
reactions, single-jet and productions in pp and pp collisions. 


5.1 Single-jet production in pp and pp collisions 

The cross section for the production of a single-jet of rapidity y and transverse 
momentum px, in a pp or pp collision is given, at lowest-order (LO), by 


d?a 1 2 /■^ XaXb 

E-r^ = > f-/ dXa- - X 

dp3 “ 1 + dij TV 2xa - Xxey 

fi{Xa,Q‘^)fj{Xb,Q^)^^{s,i,u) + (i o j) 
at 


(26) 


where xt = 2pt/\/s, Xq = xt^^/{ 2 — Xb = XaXTe~y/{2xa — xtc^) 

and ^/s is the center of mass energy of the collision. In the above sum, 
i,j stand for initial gluon-gluon, quark-gluon and quark-quark scatterings, 
ddij/dt are the corresponding partonic cross sections and is the scaling 
variable. The NLO QCD calculations were done using a code described in 
Ref. |60], based on a semi-analytical method within the ” small-cone approxi¬ 
mation, improved recently with a jet anti-fc^ algorithm for a better definition 
Ml. In Fig. im Top) our results are compared with the data from STAR 
experiment at BNL-RHIC and this prediction agrees very well with the data. 
Now we would like to test, in a pure hadronic collision, our new positive 
gluon helicity distribution, mentioned in Section 2. In a recent paper, the 
STAR experiment at BNL-RHIC has reported the observation, in single-jet 
inclusive production, of a non-vanishing positive double-helicity asymmetry 
for h <pt < 30GeV, in the near-forward rapidity region |63] . We show 
in Fig. \\M Bottom) our prediction compared with these high-statistics data 
points and the agreement is very reasonable. 


^ We thank Werner Vogelsang for providing us with the code to make this calculation. 
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Figure 14: Top: Double-differential inclusive single-jet cross section in pp 
collisions at \/s = 200GeV, versus pri^^, with jet radius parameter R=0.7, for 
—0.8 < rj < 0.8, from STAR data, obtained with an integrated luminosity of 
5.39pb“^ [62] and the prediction from the statistical approach. 

Bottom-. Our predicted double-helicity asymmetry for single-jet produc¬ 
tion at BNL-RHIC in the near-forward rapidity region, versus pt and the 
data points from STAR [63|, with the corresponding error band. 
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There are several data sets for the cross section of single-jet prodnction 
which will allow ns to test onr predictions. First we show in Fig. [T^ Top), the 
resnlts, versns for different rapidity bins, from DO [HI] and the resnlts 
from ALICE [BS] in Fig. Bottom). Except ALICE, STAR and DO are 
in good agreement with the statistical approach, as well as the resnlts from 
ATLAS and CMS displayed in Fig. [16] at y/s = 7TeV. We are giving in Table 
3 the detailed nnmber of data points and for all these experiments. Given 
the fact that the experimental results are falling off over more than eight 
orders of magnitude, this is a remarkable conhrmation of the Standard Model 
expectations, leaving no room for new physics. However some deviations 
might occur when LHC will reach a higher energy regime and this is why 
we give in Fig. HT] our predictions for the single-jet cross section at ^/s = 
13TeV. 


Experiment 

\/i TeV 


^data 

xVd-o.f. 

STAR |62] 

0.2 

11 

16 

0.67 

DO |M] 

1.96 

97 

110 

0.88 

ALICE [H3] 

2.76 

67 

20 

3.38 

ATLAS |66] 

7.0 

127 

88 

1.45 

CMS [67] 

7.0 

374 

132 

2.84 


Table 3: Detailed prediction for the inclusive single jet production. 


5.2 production in pp and pp collisions 

Let us recall that for the production in pp collision, the differential 
cross section /dy can be computed directly from the Drell-Yan picture 
in terms of the dominant quark-antiquark fusion reactions ud —)■ W~^ and 
ud —)■ W~ . So for W~^ production, we have to LO 

da^'^/dy u{xi,M^)d{x2,M^) + d{xi,M^)u{x2,M^) , (27) 

where xi^2 = Mw/ y/s exp{±y), y is the rapidity of the W and ^/s denotes 
the c.m. energy of the collision. For W~ production, we have a similar 
expression, after quark flavors interchanged and clearly these y distributions 


31 







Figure 15: Top: Double-differential inclusive single-jet cross section in pp 
collisions at ^/s = 1.96TeV, versus pP^^, with jet radius parameter R = 0.7, 
for different rapidity bins from DO |6l] and the predictions from the statistical 
approach Bottom-. Same from ALICE [65] in pp collisions at ^/s = 2.76TeV, 
with R = 0.2, 0.4 and \r]\ < 0.5. 
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Figure 16: Top: Double-differential inclusive single-jet cross section in pp 
collisions at \/s = 7TeV, versus with jet radius parameter R = 0.4, 
for different rapidity bins from ATLAS [66] and the predictions from the 
statistical approach 

Bottom-. Same from CMS [67], with R = 0.7. 
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Figure 17: Predicted cross sections for single-jet production in pp collisions 
at ^/s = 13TeV, versus p!^^, with jet radius parameter R = 0.4, for different 
rapidity bins. 


are symmetric under y —)■ —y. In the case of pp collision we have 

/dy ~ u{xi, M^)d{x 2 , M^) + d{xi, M^)u{x 2 , M^) , (28) 


which is no longer symmetric under y —>■ —y. However W~^ and W produc- 

d<j^~ 

tion cross sections are simply related since we have {y) = dy (~^)- 
Let us now turn to the charge asymmetry dehned as 


A{y) 


— PP 


dy 


dy 


iy) 




dy 


dy 


-(y) 


(29) 


and clearly we have A{y) = —A{—y). 

It contains very valuable informations on the light quarks distributions inside 
the proton and in particular on the ratio of down-to-up-quark, as noticed 
long time ago [SS]- Although the cross sections are largely modihed by NLO 
and NNLO QCD corrections, it turns out that these effects do not affect 
the LO calculation of the charge asymmetry [69]. A direct measurement of 
this asymmetry has been achieved at FNAL-Tevatron by CDF [TO] and DO 
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d] and the results are shown in Fig. [181 together with the prediction of 
the statistical approach. The agreement is very good in the low-?/ region. 
However in the high-?/ region the charge asymmetry might not flatten out, 
following the behavior of our predicted d{x)/u{x) ratio in the high-a; region 
(see Fig. 4 of Ref. [5]). 

In view of forthcoming data from the LHC, we display in Fig. [191 predictions 
from the statistical approach at ^/s = 7 and 13 TeV. In this case the W charge 
asymmetry is symmetric in yw and at hxed yw it decreases for increasing 
energy. 



Figure 18: The measured W production charge asymmetry at FNAL- 
Tevatron [701 E], versus the W rapidity yw and the prediction from the 
statistical approach, including the corresponding error band. 


However, it is not always possible to reconstruct the IT-boson and to mea¬ 
sure the boson rapidity, because of the energy carried away by the neutrinos 
in leptonic hF-boson decays. A quantity more directly accessible experimen- 
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Figure 19: The predictions from the statistical approach for the W produc¬ 
tion charge asymmetry at the LHC energies, 7TeV (solid line) and 13TeV 
(dashed line), versus the W rapidity yw 


tally is the lepton charge asymmetry, dehned as 

^ da/dr]{W+ l+u) - da/dr]{W- /~h) 

da/dr]{W+ l+u) + da/dr]{W- l-u) ’ ^ ^ 

where da/drj is the differential cross section for hF-boson production and sub¬ 
sequent leptonic decay and rj = — In [tan (6*/2)] is the charged lepton pseu¬ 
dorapidity in the laboratory frame, with 6 being the polar angle measured 
with respect to the beam axis. 

There was an earlier experimental result at the LHC from ATLAS [72], 
obtained with a total integrated luminosity of 31pb“^, but more recently 
CMS [73] released a data sample corresponding to a total integrated lumi¬ 
nosity of. 4.7fb“^. We display in Fig. [20] both data sets, together with the 
results of our calculations, which were obtained using the FEWZ code mi. 
Although the statistical approach is compatible with CMS data, a higher 
accuracy is required before considering that it is a very conclusive test of our 
PDF’s. 

Finally we consider the process —?• hF^ + X + X, where the 
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arrow denotes a longitudinally polarized proton and the outgoing have 
been produced by the leptonic decay of the hh^-boson. The parity-violating 
asymmetry is dehned as 


A 


PV 

L 


da^ — da^ 
da+ + da- 


(31) 


Here ah denotes the cross section where the initial proton has helicity h. It 
is an excellent tool for pinning down the quark helicity distributions, as hrst 
noticed in Ref. ng. 

was measured recently at RHIC-BNL [75] and the results are shown 
in Fig. 1211 As explained in Ref. [6], the W~ asymmetry is very sensitive 
to the sign and magnitude of Au, so this is another successful result of the 
statistical approach. 



Muon |? 7 | 

Figure 20: The /i charge asymmetry from IF-boson decays in bins of absolute 
muon pseudorapidity at the LHC 7TeV, with some kinematical cuts > 
20GeV for ATLAS [72] and pf^ > 25, 35GeV for CMS [73] with the predictions 
of the statistical approach (BS). 
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Figure 21: The measured parity-violating helicity asymmetries for 

charged-lepton production at BNL-RHIC from STAR [76], through produc¬ 
tion and decay of hF^-bosons versus ye the charged-lepton rapidity. The 
solid curves are the predictions of the statistical approach, including the 
error bands. 
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6 Concluding remarks 

Our quantum statistical approach to parton distributions, proposed thirteen 
years ago, has been revisited in the light of a large set of most recent world 
data from spin-averaged and spin-dependent pure DIS experiments, exclud¬ 
ing semi-inclusive DIS and hard scattering processes. The construction of 
the PDF allows us to obtain simultaneously the unpolarized distributions 
and the helicity distributions, a rather unique situation. In the current lit¬ 
erature one hnds, on the one hand, global QCD analysis of spin-averaged 
experiments, some including LHC data, [58l [771 US ES IHS E] to extract 
unpolarized distributions and, on the other hand, global QCD analysis of 
spin-dependent experiments [IS HTJ IHS IH3] , to determine the helicity distri¬ 
butions. They don’t restrict themselves to DIS experiments, like we do, and 
in general their parameterizations involve many free parameters, whose total 
number and physical meaning are considered to be irrelevant. 

Our aim in the analysis of the DIS data was to incorporate in the struc¬ 
ture of the PDF some physical principles, like Bose-Einstein distributions for 
gluons and Fermi-Dirac distributions for quarks and antiquarks, which are 
simply related from chiral properties of QCD. This allows us to reduce the 
number of free parameters, some of them having a physical interpretation. 
The improvements we have obtained from this new version is a more accu¬ 
rate determination of light quarks, strange quarks, strongly related to their 
corresponding antiquarks, and gluon distributions. We have found a large 
positive gluon helicity distribution, AG{x,Q‘^), very similar to that coming 
from the results of Ref. [82] (see also Ref. [83]) and compatible with a new 
high energy hadron collider result. As we have seen there are several chal¬ 
lenging questions related to large-x predictions, because the large-a; structure 
of hadrons is essential for a complete picture. 

The predictive power of our approach lies partly in the DIS sector, but 
mainly in the rich domain of hadronic collisions, up to LHC energies. We 
have shown that our predictions for several spin-averaged and spin-dependent 
processes at RHIC, Tevatron and LHC, are already in fair agreement with 
existing data and we expect this will be confirmed by forthcoming experi¬ 
ments. 
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